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Abstract: For 2D, symmetric, anti-parallel, collisionless magnetic reconnection, a new expression for the 
reconnection rate in the electron diffusion region is introduced. It is shown that this expression can be 
derived in just a few simple steps from a physically intuitive starting point; the derivation is given in its 
entirety and the validity of each step is confirmed. The predictions of this expression are compared to 
the results of several long-duration, open-boundary PIC reconnection simulations to demonstrate 
excellent agreement. 


I. Introduction 

At Earth's day-side magnetopause and in its magnetotail, magnetic reconnection plays a key role in 
organizing the overall structure of the magnetosphere 1, 2 . At both positions, collisionless reconnection is 
generally accepted as the appropriate model 3 . In the absence of effective collisions, it has been 
proposed that off-diagonal electron pressure tensor elements and electron inertia support the 
reconnection electric field in the electron diffusion region (EDR), thereby enabling reconnection 4 ' 6 . 
Simulation studies are available that support this explanation 7 ' 11 and models of the electron motion in 
the EDR have been constructed that explain the presence of the off-diagonal electron pressure tensor 
elements and electron inertia 12, 13 . Under the assumption that the off-diagonal tensor elements and 
inertia contribution are weak compared to the diagonal pressure tensor elements, linear approximations 
that estimate the strength of the reconnection electric field have been developed 14 ' 16 . 

In this paper, an alternate approach to an expression for the collisionless reconnection rate is presented. 
This development is limited to 2D, symmetric, anti-parallel reconnection at present. The strength of the 
electron current sheet at a reconnection site is estimated on the assumption that the average electron 
velocity there is a gyration velocity about the normal magnetic field components in the outflow jets. As 
will be shown below, this assumption leads to a notably accurate estimate of the current sheet strength. 
With this starting point, a few simple steps lead to a new expression for the reconnection rate. Assuming 
that the ion contribution to the total current sheet strength can be ignored and that the out-of-plane 
reconnection electric field strength is uniform over the EDR, an expression for the reconnection rate in 
the EDR follows 
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in which v jn is the electron inflow speed at the inflow boundary to the EDR, v eAin is the electron Alfven 

speed at the same position, L is the EDR length on the outflow axis measured in units of the local 
electron inertial length, d is the width of the EDR or of the central electron current sheet on the inflow 
axis, again measured in units of the local electron inertial length, and a = 0.45 is a dimensionless 
constant defined below. 

A long duration (500&T. 1 ) , open-boundary, driven PIC reconnection simulation has been used to guide 

the derivation of (1). This derivation and the supporting simulation results are presented below. The 
degree to which the left and right sides of (1) are in agreement demonstrates the accuracy of this 
reconnection rate expression. An additional group of related driven and undriven simulations has been 
used to more generally demonstrate the accuracy of (1). This step has shown that a is independent of 
driving rate and mass ratio, thus suggesting the possibility of a universal dimensionless constant. 

The unexpected (it/) ' dependence of the reconnection electric field that follows from (1) is in 
apparent contradiction with the behavior E rec ~ d/L that has been shown 17, 18 to follow from 

fundamental MHD conservation laws. A new expression for the EDR width d is introduced that removes 
this apparent contradiction. This expression provides the link between the conservation law approach to 
the reconnection rate and the "electron turning" approach discussed in this paper. The two approaches 
are found complementary and compatible. 


II. Simulation Setup 

A 2>2-dimensional electromagnetic PIC code described in Hesse et al. 15 has been modified to incorporate 
open inflow and outflow boundary conditions 19 ' 21 . Densities and fluxes are accumulated on the grid 
using a rectangular particle shape function. Ghost particles required by the zero-gradient boundary 
condition on the particle distributions contribute to the densities and fluxes at the open boundaries at 
each time step. The electromagnetic fields are integrated implicitly to avoid the Courant constraint on 
the propagation of light waves 15 . Lengths are normalized with respect to the ion inertial length 

d t = c co pi =c(e 2 n 0 /s 0 m i ) V using the initial current sheet density n 0 , time is normalized to the inverse 
ion cyclotron frequency co ci =eB 0 /m i using the asymptotic magnetic field strength B 0 , and velocities 
are normalized to the Alfven speed V A0 computed using n 0 and B 0 . An (x,y,z) coordinate system is 
used with the outflow and inflow directions on the x - and z-axes respectively and with out-of-plane 
directions on the v -axis. The initial equilibrium configuration is a Harris sheet 

B x = tanh(z//l) (2) 

with an additional perturbation given by 

B xp = -^-sm(^x/L x )sm(nz/L : ) (3) 
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and 


B_ p = -■^-cos(^x/Z x )cos(^z/Z_) (4) 

L x 

with A = 0.5 d i and a 0 = 0.25 , leading to a 0.5% perturbation field. No guide field is included. 

The simulations discussed below were initialized with approximately (2 -4) x 10 s particles on a 

1600x1200 grid in the xxz directions. Due to the open boundaries, these numbers vary with increasing 
simulation time. Four particle species, two of ions and two of electrons were included. The foreground 
set of ions and electrons was initialized to establish the pressure and current densities required by (2)- 
(4) and the background set was initialized to provide a constant background density n h = 0.2 . 

Background and foreground particle temperatures were initialized with T i =T e = 0.25 . Using the 
asymptotic magnetic field strength to define the electron cyclotron frequency a> ce =eB 0 /m c and the 

current sheet density n 0 for the electron plasma frequency co pe = (e 2 n 0 /s 0 m e ) l/ ~ , co pe l(o ce = 2 has been 
chosen. Further details can be found in Klimas et at. 19 ' 21 . 

The results of a group of simulations have been used to derive and test the reconnection rate expression 
(1). Of these, the primary simulation that will be used to demonstrate the derivation will be denoted 
"run60". This is a driven simulation with mass ratio mjm c = 100 on a domain L x = 80 d t by L_ = 40 d, . 
The remaining simulations were driven by varying rates, or undriven, with mass ratio m i /m e = 25 on the 
domain L x = l60d ! by L_ = 80d. . This group will be used to demonstrate the performance of the 
derived reconnection rate expression as well as to provide some evidence for the possible universality of 
the parameter a . 


III. EDR Definition 

Using the electron momentum moment equation, the EDR is defined as that region surrounding the x- 
line position in which the out-of-plane non-ideal electric field component 


E;= £ y+ [y e xB\ y 
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From (1), it can be seen that the dimensions of the EDR play an important role in defining the 
reconnection rate. These are defined in some detail in this section. 


Taken from run60 at approximately the midpoint of the run, panel (a) of Figure 1 shows E* v (z) on a cut 
through the EDR in the z -direction at the position of the x-line. The heavy vertical lines show the 
positions surrounding the x-line at which E* = 0 . These mark the boundaries of the EDR in the z- 
direction and the distance between them is the width d of the EDR. Panel (b) shows the out-of-plane 
component of the electron current j (z) on the same cut through the EDR. It can be seen that the 
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width of the EDR is also the width of the central peak in the electron current sheet. The zeroes in E* (z) 

have been used to measure the width d but, as will become evident below, the width of the current 
sheet peak is of equal importance. 



- 1.5 - 1.0 - 0.5 0.0 0.5 1.0 1.5 

Z Coordinate (shifted) 


Figure 1: (color online) (a) E* (z) and (b) j ey (z) both on a cut through the x-line position. 

The z coordinate has been shifted so that z = 0 is at the x-line position. 

Both curves have been time-averaged over a small portion of an ion gyro-period 
and smoothed with 5-point box-car averaging in z . 

The length of the EDR can be measured in the same way as the width, by finding the zeroes in E* (x) 

while moving away from the x-line in the outflow directions; the length L is defined as the distance 
between the two zeroes. A turning length also can be defined. Moving away from the x-line in an 
outflow direction, a point can be found at which the gyro-radius of a thermal electron in the outflow B. 
field is equal to the distance from that point to the x-line. This point is defined as a turning point and the 
turning length is the distance between the two turning points. Figure 2 shows the evolution of the EDR 
length over the course of run60 with 1.15 times the turning length superposed. With this multiplicative 
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factor included, the two lengths are essentially identical. Similar results have been found earlier with a 
multiplicative factor 1.3 for mass ratio 25 simulations 20 . This suggests that the two lengths may converge 


as the mass ratio is increased. For the analysis discussed below of the mass ratio 100 run60, these two 
lengths have been considered interchangeable. 


IV. Electron Turning 

This first step in obtaining the reconnection rate expression (1) is to estimate the strength of the 
electron current sheet at the x-line; i.e., the strength of the peak in j shown in Figure 1. Assuming that 

this group of electrons moves freely in the out-of-plane direction while they are turned out of the 
current sheet and into the outflow directions through gyration about the outflow B_ field, a possible 
estimate is given by 


in which j ey is the average of j ey over the width d , is the electron density at the x-line, and a) ec 
is the electron cyclotron frequency in the B_ field at a turning point (later denoted B out ). The factor 
L/2 plays the role of a gyro-radius. 
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Figure 2: (color online) Black: EDR length as defined through zeroes in E* (x) . 
Red (gray): Turning length multiplied by constant factor 1.15. 
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( 6 ) 



8 December 2014 


Page 5 



Figure 3: (color online) Black: Measured j . Red (gray): j predicted by (6). 
Both curves smoothed by running 3-point box-car averaging. 


Figure 3 shows a comparison of the measured j ey from run60 to the fit given by expression (6). It should 
be noted that there are no free parameters available to adjust the fit. 


V. Outflow Speed 

Expression (6) can be used to obtain an estimate of the electron outflow speed at a turning point; i.e., at 
the outflow boundary of the EDR. Assuming a steady state and ignoring the ion contribution to the 
current sheet strength, 
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in which B in is the magnitude of B x at an inflow boundary of the EDR and B out is the magnitude of B_ 
at an outflow boundary. Equating the second and last terms of (7) leads to 


1 n l 

B = — ^ B ,Ld = —B ,Ld 

4 m oul 4 °‘“ 


( 8 ) 


in which the carats indicate that the length and width are measured in units of the local electron inertial 
length. With v in for the magnitude of v e . at an inflow boundary to the EDR and with v mt for the 

magnitude of v at an outflow boundary, it follows that v. B =v B since £ is uniform over the 

ex ' * in in out out y 
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EDR and the non-ideal parts of E v are, by definition, zero at the boundaries of the EDR. Multiplication 
of (8) by v m leads to v in B in = v out B out = (l/4)v in B out Ld , or 


\ 

4 


Ku, =-v i „Ld 


(9) 
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Figure 4: (color online) Green (light gray): Measured v in (3-point running box-car averaged). 
Black: Measured v out (3-point running box-car averaged). 

Red (gray): v ou , fit given by (9) (9-point running box-car averaged). 


Figure 4 shows the measured v.„ and v out along with the fit to v oul given by (9). Again, there are no free 
parameters for adjustments. Notice the counter intuitive dependence on d rather than \/ d . 


VI. Reconnection Rate 

Inversion of (9) and division by v eAin leads to an expression for the reconnection rate 

v. v I 

— m_ _ 4_ou( — _ ( 10 ) 

V eAin V eAin Ld 

This expression can be simplified further. 
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Time 

Figure 5: (color online) Black: Time evolution of ratio v out /v eAin . Red (gray): Time average of v ow /v eAin 
avoiding initial transients (Time > 100). Smoothed by running 7-point box-car averaging. 



Time 

Figure 6: (color online) Black: Measured reconnection rate v in /v eAjn . Red (gray): Fit provided by 
reconnection rate expression (1). Both curves smoothed by running 7-point box-car averaging. 

Figure 5 shows the evolution with time of the ratio v out /v eAin . It can be seen that, with the exception of 
an early transient, the ratio fluctuates about an almost constant value. It should be noted that, ignoring 


8 December 2014 


Page 8 


the small difference between inflow and outflow electron densities at the respective EDR boundaries, 
the existence of this almost constant value implies 
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( 11 ) 


The electron outflow dynamic pressure at the EDR outflow boundary is approximately 1/4 the inflow 
magnetic field pressure at the inflow boundary. Below, evidence that suggests this condition is generally 
true will be discussed. 


Substitution of 


a = 


_ time averaged 


( 12 ) 


for the ratio v out /v eAin in (10) leads to the reconnection rate expression (1). Figure 6 shows a comparison 
of the right to the left sides of (1) evaluated over the course of the run60 simulation with a = 0.45 , the 
time-averaged value with the initial transient evident in Figure 5 excluded from the average. Again, it 
should be noted that there are no available free parameters for adjustments of this fit. The source of the 
overshoots in the fits to the peaks in the measured reconnection rate is under investigation. 


VII. Out-of-Sample Performance 

The parameter a plays a prominent role in normalizing the reconnection rate expression (1). If a were 
to depend sensitively on specific details of a reconnection scenario, then the utility of (1) might be 
questioned. Some evidence suggesting a universal value for a has been found. 


Simulation 

m i / m e 

Driver Strength 

Run Time yco cl . j 

a 

run60 

100 

0.3 

500 

0.45 

run59 

25 

0.3 

500 

0.47 

run58 

25 

undriven 

372 

0.49 

run54 

25 

0.4 

520 

0.46 

run53d 

25 

0.2 

500 

0.46 


Table 1: (color online) a values for several simulations of various driving rates and two mass ratios 

Table 1 shows the values for a that have been found for a group of simulations, including run60. In all 
cases a was obtained through the time average (12) with initial transient behavior excluded. Driver 
Strength is the value of the out-of-plane electric field imposed at the inflow boundaries of a simulation 
following a short start-up interval during which the field strength was ramped up from zero. Comparison 
of run60 and run59 shows the effect of varying the mass ratio while holding the driver strength fixed. 
Simulations other than run60 show the effect of varying driver strength holding the mass ratio at 25. For 
the undriven run58, the out-of-plane field at the inflow boundaries was allowed to float with a zero 
gradient imposed. Admittedly, Table 1 contains a small sample but it is difficult to ignore the suggestion 
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that a represents a universal relationship between electron outflow dynamic pressure and inflow 
magnetic field pressure as expressed by (11). 




100 200 300 400 500 100 200 300 400 500 

Time Time 

Figure 7: (color online) Black: Measured reconnection rate. Red (gray): Fit given by (1) 
for (a) run59, (b) run58, (c) run54, and (d) run53d using respective values for a from Table 1. 

All curves smoothed using running 3-point box-car averaging. 

The performance of the reconnection rate expression (1) for run60 is shown in Figure 6. Figure 7 
continues with examples of that performance for the remaining mass ratio 25 runs listed in Table 1. 
Again, there are no free parameters available and, again, there are occasional problems with excessive 
spikes in the fits. The source of these spikes is under investigation. 


The degree to which the fits shown in Figure 7 succeed in reproducing the measured reconnection rates 
was unexpected. Two assumptions that were made in the derivation of (1) do not strictly apply. First, a 
cut in the z variable through a mass ratio 25 EDR, as in Figure 1, generally yields a shape that is similar 
to that of the current distribution shown in panel (b) of Figure 1 (e.g., see Figure 6, Klimas et a I. 20 ). A 
central peak in the EDR, whose width coincides with that of the central peak in the current distribution, 
sits on top of a broader distribution that extends out to ion scales. As in the analysis of run60, the widths 
of these central peaks have been measured to fix the width d . Consequently, the non-ideal part of the 
out-of-plane electric field is not zero at the inflow boundary of the EDR, so defined, and the assumption 
leading to (9) that E at the inflow boundary is given by v jn B in fails somewhat. Second, neglect of the 
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ion current contribution to the strength of B in in (7) is less justified in the lower mass ratio simulations. 
The ion current contribution, averaged over the EDR width, can amount to as much as 15-20 percent of 
the electron contribution for the mass ratio 25 simulations; the ion contribution for the mass ration 100 
run60 is typically around 5 percent. Despite these faults, the simple expression (1) yields surprisingly 
accurate fits to the measured reconnection rates for the mass ratio 25 as well as the mass ratio 100 
simulations. 

VIII. EDR Width 

A straight forward derivation, based on fundamental conservation laws, shows that the reconnection 
electric field E rec ~ d/L 17 ' 18 . A simplified portion of this derivation is given in this section. From (10), 
however, it also will be shown in this section that E rec ~ l/Ld . Resolution of this apparent contradiction 

imposes a new relationship that restricts the value of d . The following applies strictly to run60 and 
approximately to the rest of the simulations listed in Table 1 for the reasons discussed in the preceding 
paragraph. 

Starting with a box surrounding the EDR of length L and width d , the length and width of the EDR, 
particle conservation for the electrons leads to 

{ v ,„"i„) L = { v o l ,t n o.«) d (!3) 

in which the brackets ( ) indicate averages over the respective box edges. A detailed examination has 
shown that {v in n in ) = (v in )(n in ) = v in ri in quite accurately. In addition, {v out n out )^{v out ){n out )^{v out )n out but 
( v out)^ v out ■ < 3 )n the outflow edges of the box, the outflow speed is peaked in the centers of the outflow 
jets such that v out = 2(v out } . From (13) 


v B. =(v ,)B. - = E (14) 

in in \ out / in rec ' 

in which E is the reconnection field strength and n ,ln. - 1 has been assumed. 

rec out / in 

From (10) it can be seen that 


I TYl 1 
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(15) 


Equating (14) and (15) leads to 



(16) 
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Time 


Figure 8: (color online) Black: Measured normalized EDR width. Red (gray): Fit given by (16). 

Both curves smoothed using running 7-point box-car averaging. 

Figure 8 shows that the normalized EDR width does satisfy (16). Thus, the independent expressions (14) 
and (15) for the reconnection electric field strength are complementary and compatible, as are the 
independent methods used to obtain those expressions. Similar results have been obtained for all of the 
simulations listed in Table 1. Although, as in Figure 8, each of the simulations shows some variability in 
time, generally it has been found that d- 3, independent of the two mass ratios and various driving 
rates of the group. 


IX. Discussion 

For symmetric, anti-parallel, collisionless reconnection a new expression for the reconnection rate has 
been derived in just a few simple steps from a physically intuitive starting point, the expression (6) for 
the electron current sheet strength. That expression shows that the average motion of a free electron 
at the x-line position in the current sheet should be viewed as the initial phase of a gyration about B_ at 
a turning point. Surprisingly, the strength of the reconnection electric field does not play a direct role in 
fixing the current strength in (6). 

In addition to the expression (6) for the electron current sheet strength, the derivation of the 
reconnection rate expression yields two interesting new results before the final expression is reached. 
The first of these is the expression (9) for the electron outflow speed at the EDR outflow boundary. A 
commonly employed concept of plasma conservation for plasma entering the EDR through the inflow 
boundary and then leaving through the outflow boundary suggests that, holding other parameters fixed, 
the outflow speed should depend inversely on the width d ; as the width of the outflow boundary 
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decreases outflow speed must increase to maintain a steady throughput. However, (9) states that the 
outflow speed is proportional to d and Figure 4 confirms the accuracy of this expression for the outflow 
speed. This apparent contradiction, as well as the apparent contradictions indicated by the expressions 
(14) and (15) for the reconnection electric field strength, is removed by the restriction on d given by 
(16). For the tightly coupled nonlinear reconnection system, once the system is specified, then the width 
d is no longer free for adjustment. Plasma conservation applies but the commonly employed 
adjustability of d does not. 

The second new feature discovered on the way to the reconnection rate expression is the quasi-steady 
evolution over almost the entire run60 of the ratio v ou Jv eAin as exhibited in Figure 5, which leads to the 
approximate relationship between outgoing electron dynamic pressure and inflowing magnetic pressure 
given by (11). The final reconnection rate expression (1) is obtained by replacing this ratio by its time 
average a . As exhibited in Table 1, a appears to be independent of mass ratio as well as driving rate 
for driven and undriven simulations. This unexpected result further increases the significance of (11) and 
also suggests an investigation of the invariance of a under additional changes in the reconnection 
scenario. 

If the universality of a can be demonstrated, then all of the terms in the reconnection rate expression 
(1) except for the EDR length L should be measurable in Earth's magnetotail by NASA's upcoming 
Magnetospheric Multiscale (MMS) mission (http://mms.space.swri.edu/index.html). Instead of I, it is 
possible that the equivalent turning length may be accessible from MMS measurements if three or four 
of the spacecraft pass through a single EDR sufficiently closely in time so that the dependence of the 
outflow field strength on displacement in the outflow directions can be modeled. If so, then the 
positions of the turning points may be modeled and the turning length estimated. Given that it will be 
possible, in this multi-spacecraft mission, to maintain the inter-spacecraft spacing at as small as 10 km 
the possibility of multiple spacecraft EDR passages is not out of the question. 

In addition to observational verification of (1), if it is possible, further reconnection simulation studies 
are necessary. All of the results presented in this paper are noisier than one would wish; simulations 
with higher macro-particle densities should be carried out. The suggested convergence of the EDR and 
turning lengths with increasing mass ratio should be further tested. Verification would increase the 
likelihood of observational confirmation as discussed in the preceding paragraph. Confirmation of all of 
the results presented in this paper requires simulations using other available codes. Based on this 
experience, however, it is firmly believed that true open-boundary codes will be necessary. 

Finally, the simplicity of the derivation discussed in this paper suggests that further extensions to the 
more realistic anti-symmetric and/or guide field reconnection models may be possible. 

Acknowledgments: This research was supported by NASA's MMS IDS grant NCC5-494 (MOST). 

Invaluable discussions with Paul Cassak are greatly appreciated. 


8 December 2014 


Page 13 


X. References 

X A. Nishida, Space Sci. Rev. 91 (3-4), 507-577 (2000). 

2 V. M. Vasyliunas, Space Sci. Rev. 158 (1), 91-118 (2011). 

3 R. A. Treumann, W. Baumjohann, arXiv:1401.5995vl [physics. space-ph] (2014). 

4 V. M. Vasyliunas, Rev. Geophys. 13 (1), 303-336 (1975). 

5 J. W. Dungey, La Recherche Spatiale 2 (15) (1989). 

6 M. Hesse, T. Neukirch, K. Schindler, M. Kuznetsova and S. Zenitani, Space Sci. Rev. 160 (1-4), 3-23 
( 2011 ). 

7 M. Hesse and D. Winske, J. Geophys. Res. 99 (A6), 11177-11192 (1994). 

8 R. Horiuchi and T. Sato, Phys. Plasmas 4 (2), 277-289 (1997). 

9 M. M. Kuznetsova, M. Hesse and D. Winske, J. Geophys. Res. 105 (A4), 7601-7616 (2000). 

10 J. Birn, J. F. Drake, M. A. Shay, B. N. Rogers, R. E. Denton, M. Hesse, M. Kuznetsova, Z. W. Ma, A. 

Bhattacharjee, A. Otto and P. L. Pritchett, J. Geophys. Res. 106 (A3), 3715-3719 (2001). 
n H. Karimabadi, W. Daughton and J. Scudder, Geophys. Res. Lett. 34, L13104 (2007). 

12 K. Fujimoto and R. D. Sydora, Phys. Plasmas 16 (11) (2009). 

13 A. Divin, S. Markidis, G. Lapenta, V. S. Semenov, N. V. Erkaev and H. K. Biernat, Phys. Plasmas 17 (12) 
( 2010 ). 

14 M. M. Kuznetsova, M. Hesse and D. Winske, J. Geophys. Res. 103 (Al), 199-213 (1998). 

15 M. Hesse, K. Schindler, J. Birn and M. Kuznetsova, Phys. Plasmas 6 (5), 1781-1795 (1999). 

16 M. M. Kuznetsova, M. Hesse and D. Winske, J. Geophys. Res. 106 (A3), 3799-3810 (2001). 

17 P. A. Cassak and M. A. Shay, Phys. Plasmas 14 (10) (2007). 

18 P. A. Cassak and M. A. Shay, Phys. Plasmas 16 (5) (2009). 

19 A. Klimas, M. Hesse and S. Zenitani, Phys. Plasmas 15, 082102 (2008). 

20 A. Klimas, M. Hesse, S. Zenitani and M. Kuznetsova, Phys. Plasmas 17, 112904 (2010). 

21 A. Klimas, M. Hesse and S. Zenitani, Phys. Plasmas 19, 042901 (2012). 


8 December 2014 


Page 14 


